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ABSTRACT 

This paper deals with Automatic generation control of two area interconnected hydro -thermal system using 
proportional plus integral controllers with Particles swarm optimization techniques (PSO). The parameters of 
PI controllers are tuned with the help of optimal control strategies such as Integral square error (ISE) and 
Integral time multiplied hy absolute error (ITAE) by particles swarm optimization technique. Mismatch in the 
total generation to total demand plus losses results in the system frequency deviation from its nominal value and 
causes variation in the Tie-line power flow. System performance is examined considering perturbation of .1 pu 
step load change in thermal system. Simulation studies and a comparison of dynamic responses reveal the 
performance of particle swarm optimization in terms of the frequency and Tie-line power responses in both the 
areas. Comparison is also made between the two different performance indices criteria (ISE, ITAE). 

Keywords - Automatic Generation control, Particle swarm optimization, PI controllers and Performance 
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I. INTRODUCTION 

An interconnected power system consists of different areas which are connected to each other by Tie -lines. In 
power system both the active and reactive demands are never the same and they continuously changes with time. 
Automatic Generation Control (AGC) of interconnected power systems is gaining importance in modern power 
systems. Mismatch between the total generation to total demand and associated system losses, results in 
deviation of system frequency from its nominal value and also changes in the Tie-line power flow from the 
desired value. In each area all the generators changes their speed to maintain the frequency and change in Tie - 
line power to the desired value in static as well as in dynamic conditions. The primary objectives of a control 
area are to satisfy its local consumer demands and to respond to the variation in demand of other control areas. 
Frequency and tie-line power exchange are the two variables of interest, and combination of these two called 
area control error (ACE). Automatic generation control (AGC) is the name given to a control system having 
three major objectives [1,2, and 3]: 

a. To hold system frequency at or very close to a specified nominal value. 

b. To maintain the correct value of interchange power between control areas. 

c. To maintain each units generation at the most economical value. 

A literature survey shows that lot of research paper has been carried by linearized model of thermal as well as 
hydro units in both the single area and in interconnected areas [4, 5]. Most of the works concerned with AGC of 
interconnected power systems pertain to tie-line bias control strategy [6]. In past tuning of parameters of 
different controllers (PID, PI) is done by various methods like maximum peak resonance specification (MPRS) 
[7]. AGC problem for hydro thermal systems have also studied in [8,9, 10]. Tuning the parameters of PI 
controllers in the two area hydro-thermal system connected to each other are done by Genetic algorithm and 
other optimization techniques like PSO, and fuzzy logic control [11-12]. 

The objective function which is minimized by the use of optimization technique is the different performance 
indices criteria. These criteria are Integral of the squared error (ISE), Integral time multiplied by the absolute 
error (ITAE). 

oo 

IS E = J e 2 (t)dt (1) 


ITAE = j t \ e (O I dt 
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In this paper Particles swarm optimization (PSO) is used to tune the parameters of the PI controllers in the 
hydro-thermal system. The comparison is made between the two performance indices criteria of dynamic 
responses of the deviation in frequency and change in Tie-line power flow in terms of undershoot and settling 
time. The parameters which are tuned by optimization techniques (PSO) are Kprl, Kil, Kpr2 and Ki2. This 
paper is described as follows: 

PI controller described in section II. In the section III State space modelling is carried out. In section IV test 
system is explained. Section V includes the explanation and application of PSO. Results and analysis are given 
in section VI. Section VII concludes the paper 

II. PIDCONROLLER 

Output of PID controller is given by: u ( t ) = k pr e(t)+ K.\ e(t)dt (3) 

Where e ( r ) is error input to PI controller i.e. area control error, u(t) is the output of PI controller, K 
proportional gain and K . integral gain. 


III. STATE SPACE MODELING 

3.1 Transfer Function 

A two area hydro-thermal interconnected power system is considered to tune the parameters of the PI 
controllers. Thermal area used is non-reheat turbine unit connected with hydro unit with supplementary 
controllers. Perturbation of 0.1 pu step load change is considered in thermal area. Transfer function model of the 
two area hydro-thermal system is shown in Fig (1). PI controller is used for the automatic generation control of 
the system. These PI controllers are used to bring the system frequency to the nominal value and also the Tie- 
line power to the desired value. The transfer function of the desired relation can be obtained by state space 
modeling of the system. Two area hydro -thermal systems can be described by standard state space equations 


X = /\ X -b B ZJ 



y = cx + n u (5) 

Where U is the input vector, X is the state vector and Y is the output vector. And A, B, C, D are the real constant 
matrices of appropriate dimensions which in turn depend on the system parameters. These vectors in equation 
(4) and (5) can be represented as: 
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State space matrix can be obtained by the use of differential equations. State equations for the different 
components of the interconnected power system can be shown by the following equation: 
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By the use of these differential equations state space matrix can be find out. These real constant matrices are: 
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With the help of these matrices, we get transfer function of the system as defined: 

Transfer function = C[S/-A] B + D 

This transfer function is in the matrix form of order (3 x 2). 
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3.2 Objective function 

For tuning of the parameters of the PI controller’s different objective functions are used and then comparison is 
made so that best system response can be obtained. These different objective functions are: 

00 
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IV. TEST SYSTEM 

The two area hydro-thermal interconnected is considered with PI controllers. One area is non-reheat thermal 
system and other area is hydro system. Both the areas are having generating equipment with different parameters 
[13]. Simulation is carried out with the step load change of .1 pu in thermal area. This will affect the frequency 
response and tie-line power change. For controlling this PI controllers are used. Transfer function model of two 
area hydro-thermal interconnected system is shown in Fig. 1 . 



Fig 1 . Transfer function model of hydro-thermal interconnected system 


V. PARTICLES SWARM OPTIMIZATION 

Particle swarm optimization (PSO) is a population based stochastic optimization technique developed by Dr. 
Eberhart and Dr. Kennedy in 1995, inspired by social behaviour of bird flocking or fish schooling [15]. The 
system is initialized with a population of random solutions and searches for optima by updating generations. 
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However, unlike GA, PSO has no evolution operators such as crossover and mutation. In PSO, the potential 
solutions, called particles, fly through the problem space by following the current optimum particles. Each 
particle keeps track of its coordinates in the problem space which are associated with the best solution (fitness) it 
has achieved so far. (The fitness value is also stored.) This value is called pbest. Another "best" value that is 
tracked by the particle swarm optimizer is the best value, obtained so far by any particle in the neighbors of the 
particle. This location is called lbest. When a particle takes all the population as its topological neighbors, the 
best value is a global best and is called gbest. A binary PSO algorithm has been developed in [16]. In the binary 
version the particle position is not a real value, but either the binary 0 or 1 . The velocity is used as a probability 

threshold to determine whether x ^ , the dth component of x * should be evaluated as a zero or a one. They 

squashed x id in a logistic function 


1 

s ( x id > = 

1 + exp(-v„) 


(24) 


Then generated a random number r for each bit-string site and compared it to s(x . d ) . If r was less that the 
threshold, then x was interpreted as 1, otherwise as 0. The velocity of the each particle can be updated by the 

id 

equation given below: 

v./ +1 = v./ + rl* cl(pbest t - x.f ) + r2 * c2(gbest T - x.f) 

id id v r ,d 7 v o id 7 

pbest : Personal best of particle i 
gbest: Global best (best of pbest of the group) 
rl, r2: Random number in the [0, 1] interval 
cl, c2: Acceleration coefficient 

x id Position of i th particle in dimension d 


v id Component in dimension d of the i th particle velocity in iteration t 

Particle velocity is calculated by the formula shown previously and then with the help of random number, 
particles position is updated. 

Binary Particle swarm optimization is used to minimize the objective function. . Procedure used for tuning of PI 
controllers with the help of BPS O is shown below: 

1. Generate the initial particles position for each variable Kprl, Kpr2, Kil and Ki2. Here 24 population size is 
used. 

2. Find the value of the objective function (equation 22-23) with the help of each particles position, do it for all 
population and then find pbest and gbest. 

3. Calculate the velocity of each particle with the help of equation (25). 

4. Again find each particles position by equation (24). 

5. Again evaluate the objective function with the help of new particles position and repeat step2, 3 and 4. 

If no. of iteration > max generation then stop, and best solution formed, else go to next iteration. 


VI. RESULT AND DISCUSSIONS 

Simulation has been done by taking the system given in Fig. 1 with one thermal area and other hydro area. 
Particle swarm optimization techniques used to tune the parameters of the PI controllers of both the area when 
there is .1 pu load perturbation occurs in area-1. The optimal criterion is taken as ISE (integral squared error) 
and ITAE (integral time multiplied by absolute error). The simulation results for the change in frequency in both 
areas and change in the Tie-line power flow for all performance indices criteria with Particles swarm 
optimization is properly analysed. Dynamic responses of interconnected hydro-thermal power system for ISE 
and ITAE are shown in Fig. 3, 4. 

The dynamic responses of , A*: and for step load increase in thermal area with BPSO is shown in Fig. 
3 and Fig. 4, and with Table 1 it is clear that ITAE shows good dynamic response than ISE in terms of settling 
time. 

Table shows the parameters of PI controllers tuned by the BPSO. From the table it is clear that the objective 
function value ITAE is less than ISE. This is because ITAE is more selective performance index than other 
indices. Even though ITAE does not penalize large initial errors that are unavoidable but it does penalize large 
duration transient [17]. 
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Table 1 tuning of PI controllers from BPSO with ISE and ITAE 




Kprl 

Kil 

Kpr2 

Ki2 

PI 

iPLl 

Settling 

time(s)fl 

Undershoot 

Binary Particle 

ISE 

.470 

.941 

0.00 

0.0 

17.54 

.1 

60 

.2 

swarm 


6 

2 


0 

8 

p.u. 



optimization 

ITAE 

.117 

.705 

.941 

0.0 

1523. 

.1 

30 

.22 



6 

9 

4 

0 

0 

p.u. 




Table shows the parameters of PI controllers tuned by the BPSO. From the table it is clear that the objective 
function value ITAE is less than ISE. This is because ITAE is more selective performance index than other 
indices. Even though ITAE does not penalize large initial errors that are unavoidable but it does penalize large 
duration transient [17]. 



Fig 2. Flow chart of BPSO 
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Fig. 3 Dynamic responses of frequency and power by BPSO with ISE when .1 pu step load change in area-1 
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Fig.4 Dynamic responses of frequency and power by BPSO with ITAE when .1 pu step load change in area-1 


VII. CONCLUSION 

Load frequency regulation characteristic of two different hydro -thermal area power system have been studied. 
Particle swarm optimization (PSO) has been successfully applied to tune the parameters of the PI controllers. 
Two performance indices criteria are taken (integral square error (ISE), integral time multiplied by absolute 
error (ITAE). Binary Particles swarm optimization (BPSO) has been successfully applied to tune the parameters 
of the PI controllers for hydro-thermal system. From the above results it concluded that for automatic generation 
control of hydro-thermal system ITAE performance indices criteria gives better dynamic performance in 
comparison with ISE for PSO controllers. Damping responses of ITAE are good in terms of under-shoot and 
settling times than ISE. 
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VIII. Appendix 

t j =10 


Typical values for the system parameters are [13]: 
K , = K =120 Hz/p.u. MW 

pi pi ' ~ 

T = T = 20 sec 

pi p 2 

R { = R 2 = 2.4 Hz/p.u. MW 
= B 2 = .4249 

T =.08 sec 

8 

T t =.3 sec 
T l2 =0.0866 
T x = 48.709 
T 2 =.51308 


T =1 

w 

P = P , =1200 MW 

r 1 r 2 

D { = D 2 =.00833 p.u. MW/Hz 

Parameters of BPSO used to tune the PI controllers 
are shows below: 

Acceleration coefficient cl=1.5; c2=2.5; 

Damping coefficient =.99; 

No. of iteration=100; 

No. of Population=24; 
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